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Abstract Baryonic physical processes could leave non-negligible imprint on cosmic mat- 
ter distribution. Series of high resolution simulation data sets with identical initial condition 
are employed for count-in-cell (CIC) analysis, including one N-body pure dark matter run, 
one with adiabatic gas only and one with dissipative processes. Variances and higher order 
cumulants S n of dark matter and gas are estimated. It is found that baryon physical pro- 
cesses mainly affected dark matter distribution at scales less than l/i _1 Mpc. In comparison 
with the pure dark matter run, adiabatic process alone strengthens variance of dark matter by 
~ 10% at scale 0.1/i _1 Mpc, while 5„s of dark matter deviate only mildly by a few percent- 
ages. Dissipative gas run does not differ much to the adiabatic run in dark matter variance, 
but renders significantly different S n parameters of dark matter, bringing about more than 
10% enhancement to S3 at 0.1/i _1 Mpc and z = and being even larger at higher redshift. 
Distribution patterns of gas in two hydrodynamical simulations are prominently different. 
Variance of gas at z = decreases by ~ 30% in adiabatic simulation while by ~ 60% in 
non-adiabatic simulation at 0.1/i _1 Mpc, the attenuation is weaker at larger scales but still 
obvious at ~ 10ft- _1 Mpc. S n parameters of gas are biased upward at scales <^ 4/i _1 Mpc, 
dissipative processes give ~ 84% promotion at z = to S3 at 0.1/i _1 Mpc in contrast with 
the ~ 7% change in adiabatic run. The clustering segregation between gas and dark matter 
could have intricate implication on modeling galaxy distribution and relevant cosmological 
application demanding fine details of matter distribution in strongly nonlinear regime. 

Key words: cosmology: dark matter — large-scale structure of universe — methods: statis- 
tical 



1 INTRODUCTION 

The present clustering pattern of large scale structures on cosmological scales is generally interpreted as the 
growth of primordial density fluctuations mainly through gravitation instability of the dark matter which 
dominates the matter content of the Universe. While at large scales the gravitation monodrama of dark mat- 
ter is much appreciated, at small scales ignoring non-gravitational effects associated with galaxy formation 
would induce considerable systematics to relevant application. Quantification of such impact of baryons, 
including acting scale ran ge and strength, is strongly desired to mee t the accuracy budget of cosmological 
parameter estimation (e.g S haw et alj, 2010HSemboloni et all|20TTI) and structure formation and evolution 

etal.Ll2009h . 



model refinement (e.g. lStanek et al.l 2009 : Dolag et al. , 20091) . In contrast to the simplicity of gravitational 



force, physical processes baryons involved in, such as radiative cooling, and star formation etc., are usually 
very complicated and highly entangled, even worse is that baryon physics works normally in strongly non- 
linear regime where gravitational evolution is already analytically intractable. Advanced computational fa- 
cilities and algorithms, in together with accumulated knowledge summarized from modern observation, has 
enabled high resolution hydrodynamic simulations with various treatment prescriptions for baryon physics 
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plugged in (e. g. lTeyssieil l2002tlSpringelI 120051) . So to date investigati on on effects of differen t baryon phys- 
ical processes is carried on mainly with numerical simulations (e.g. Ivan Daal en et all 1201 lb though there 
are still a long way to go to build a trustworthy machinery to capture the messy gas physics in full details. 

Recent analysis of simulation data sets has shown that baryonic physical processes could alter the matter 
power spectrum at k > l OfoMpc" 1 and th en the weak lensing power spectrum consequently by some 
non-negligible percentage dJing et all 120061), which is confirmed and much extended in later works (e.g. 



Rudd ^tajiHooilHearin & Zentner, 2009; van Daalen et alll201ltlSemboloni et alll20TTtlCasarini et all 



2012). The modulation to power spectrum at such small scales, speaking in terminology of halo mode for 
matter clustering, is mainly happened upon the on e-halo term which is de termined just by mass distribution 
inside halos and halo mass distribution function dCoorav & Shethl |2002b. P r esenc e of gas in simulations 
does slightly boost concentration parameter of halo mass pr ofile (Lin et al., 2006) and regulate the high 
mass branch of halo mass function measurably dStanek et all I2009T Cui et al. , 201 lb . even with adiabatic 



process alone. S uch changes in halo prope rties actually can bring about significant amendment to strong 
lensing statistic s (Wambsganss et al., 2008) and the thermal and kinetic Sunyaev-Zel'dovich power spectrua 
dBattaglia et alll2010tlShaw et alll2010b as well. 

The basic scenario here is that baryons are directly redistributed by adiabatic contraction, radiative 
cooling, various feedbacks from galaxies and their central black holes, and star formation activities etc., 
and then dark matter distribution is modified through the gravity coupling with baryons. Since dark matter 
and baryons experience different actions, it is normal to expect that their clu stering would dif fer with the 
pattern shown in the dark matter only case in complex way. In the work of Ijing et al.l d2006l) it is found 
that the clustering of the gas is suppressed while that of dark matter is boosted at scales k > l/iMpc -1 , 
resulting in to the clustering of total matter suppress at the level 1% at 1 < k < 10/iMpc -1 but then boost 
up to 2% in the nonradiative r un and 10% in the run with star formation at k ps 20ft.Mpc~ 1 . The extensive 
research by Ivan Daalen et all (1201 ll) with AGN feedback provided a quantitative different descri ption, due 
to thei r different cooking recipes for gas physics, though it is still qualitatively in agreement with Jing et al. 
d2006h . They discovered that the 1% level decreasement of power spectrum of total matter at z — to that 
of pure dark matter simulation can be as low as k ~ 0.3/iMpc" 1 , 10% dropping appears at k ~ 10/iMpc^ 1 , 
clustering enhancement is observed somehow at k >^ 70/iMpc -1 . Compared with the power spectrum of 
pure dark matter run, gas in hydrodynamic simulations exhibits much less power for k > 1/iMpc -1 while 
dark matter component shows power boost at k > 10/iMpc^ 1 . 

The baryon influence on matter clustering and the baryon-dark matter segregation in clustering can 
be better inspected with higher order correlation functions which is known to be able to reveal more sub- 
tle details of clustering than two-point statistics like the power spectrum. iGuillet et al.l (120101) measured 
skewness of the MareNostrum simulatiorQ and also a set of dark matter only simulation. They found that 
relative to the matter distribution in pure dark matter simulation, in hydrodynamic simulation the dark mat- 
ter component has skewness decreased mildly between 0.3 < r < l/i _1 Mpc and then boosted apparently 
at smaller scales. They demonstrated that by adding an exponential gaseous disk profile to the halo model 
could roughly reproduce their measurements. 

In this paper we perform count-in-cell measurements of N-body/SPH simulations in together with a 
pure dark matter simulation as reference, in order to better depict the impact of gas physics on matter 
distribution at higher orders in compl ement to works bas ed on power spectrum, and which also serves as 
an independent check to the results of IGuillet et al.l d2010b . Furthermore, our interests is particularly on the 
gas-dark matter segregation phenomenon, i.e. distribution differences among matter of different species, 
which might cast light on the origin of the galaxy biasing. The layout of the paper is as following, section 
2 contains description to the simulation data and estimation method for higher order correlation function, 
results and their analysis are in Section 3, the last section is of summary and discussion. 



http://astro.ft.unam.es/~marenostrum 
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2 COUNT-IN-CELL MEASUREMENTS OF SIMULATION DATA 
2.1 The Simulation data sets 

The simulations data sets we use are the same as in IJing et all (120061) . which consist of three simulations 
produced by the GADGET2 code ( Springel, 2005), one pure dark simulation, one hydrodynamic simulation 
with adiabatic process alone, and one hydrodynamic simulation incorporated with radiative cooling, star for- 
mation, and supernovae feedback etc. . The three simulations are run with 512 3 particles for each component 
of dark matter and gas, starting at Zi n i = 120 with the same initial condition in a cubic box of 100/i _1 Mpc, 
and their cosmology parameters are set t o (n m , Oa , f2h, ct k, n, / t)=(0.268 , 0.732 , 0.044, 0.85, 1, 0.71. More 
details of simulations can be found in Ijing et alj d2006) and iLin et al.l d2006l) . Here three snapshots at 
z = 0, 0.526, 1.442 are picked up for analysis. 



2.2 Count-in-cell and higher order correlation function 

Here higher order correlation functions are those volume averaged ones, i.e. higher order connected mo- 
ments of smoothed density fluctuation fields 5 by certain window function w, £ n = (S n ) c - can be 
estimated through the count probability distribution function P/v (R) by count-in-cell method. Given cubic 
cell of side size R, Pm at this scale is the probability that a randomly thrown cell in the catalog contains N 
galaxies, 

1 C 

Pn = cT, Sd ^ = n ^- W 

i—l 

Under the usual local Poisson approximation, Pn is the probability d istribution function jp(<5) of smoothed 
density fluctuation convolved with a Poisson kernel (see the review of iBernardeau et ail 120021) 

where { N) is the mean count-in-cell. Higher order correlation functions £ n are often expressed by the higher 
order cumulants hierarchy S n , 

Sn = =|=T ! (3) 
S2 

which can be derived from the following recursion relation 

„ _ l 2 F n _ 1 f n \ (n - k)S n - k F k 
c k=l v 7 c 

where N c — (N)^ an d the factorial moments 

F k = J2 Pn ^ x ( N h = ( N ( N -l)---(N-k + l)) . (5) 
Explicit estimators for the variance, skewness and kurtosis are just 

- _ F 2 

e 2 - ps - 1 



s, 
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F 1 {F 3 _3F 1 F 1 + 2F?) 
F?(F 4 - 4F 3 F 1 - 3F| + 12F 2 F? ~ 6Ff) 



S 3 - 77? ( 6 ) 



{F2-m z 

Pn is calculated with the over-sampling algorithm to reach a ~ 10 7 sampling rate (Szapudi, 1998). 
Probing scale R is limited within (0.1, 10)/i -1 Mpc, the small scale cut is chosen so to ensure robust re- 
covery of statistics over discreteness (normally (N) > 0.1 is sufficient), and large scale limit comes from 
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Fig. 1 £ 2 and S n — £, n /£,2 up to n = 6 of dark matter in simulation runs of pure dark matter 
(left column), with adiabatic gas physics only (middle column), with star formation and other gas 
physics (right column). Solid, dotted and dashed lines are of z = 0, 0.526, 1.442 respectively. 



one tenth of the box size above which correlation functions are no longer reliable by practical experiences. 
Since all simulations are evolved from the same initial condition, in the same volume and of the same reso- 
lution, there is no need to calculate their cosmic variance (or error bars) if we are simply interested in their 
differences. 



3 INFLUENCE OF BARYONIC PHYSICS ON CLUSTERING 
3.1 clustering of dark matter 

Correlation functions of dark matter in three simulation runs are illustrated in Figure [1] showing analogous 
redshift evolution history. It has been checked that if we rescale £ 2 with D(z) the growth rate of large 
scale structure (see Lah av et all 1 1 99 U for approximate formula), variances at different redshift are in good 
agreement at scales where £ 2 (R> Z = 0) < 1, at smaller scales £ 2 (z)D 2 (z — 0)/D 2 (z) becomes lower 
with increasing z as what is well known. 

S n parameters are apparently larger at higher redshift at scales R <~ 3ft. _1 Mpc, the decrement from 
z = 1.442 to z = 0.526 in S n parameters is much bigger than that from z — 0.526 to z — 0. Recall that 
£ n>2 s are measures of non-Gaussianity and our simulations are evolved from Gaussian initial condition, 
intuitively non-Gaussianity would increase when redshift decreases, in line with the growth of gravitational 
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nonlinearity. Let the cosmic scale factor a = 1/(1 + z) and Act > be a small increment to ct, as S n (a + 
Act) < S n (a), S n — > £«( a + Act) > £«( a ) an d at small scales £ n > 0, there is 

logUa + Aa) - logUa) < (n - 1) [bgl a (a + Aa) - log? 2 (a)] , (7) 

subsequently 

— <(n-l)— , fom>2, (8) 

a log a a log a 

which establishes an interesting relation between the evolution rates of higher order correlation functions 
and that of the two point correlation function in strongly nonlinear regime. 

At scales R >^ 3/i _1 Mpc 5V,, s show little redshift depe ndence, there is S n (a < 1) « S n (a = 1) 
as perturbation theory predicted (Eos alba & Gaztafiagai [?998 ). However the approximation is not perfect 
in our raw results, at large scales there are differences at different epochs at level of several percentages, 
S n s at higher redshift turns to be slightly larger. It is known that there is blemish rooted in the Zel'dovich 
approximation based initial condition generator, the resulting correction to measured S n s at leading order 
decays with redshift at rate roughly oc \D{z) / Djzjnj)] -1 . But the effect of the bias is that true 5„s are 
larger than measured values (IScoccima rro, 1998; Fos alba & Gaztanagai 1 1998b . so that the offsets between 
different redshift will be even higher than what is shown in FigureQ] The puzzle could be dynamical rather 
than systematical biases, here we just leave the issue to future work. 

Influence of gas physics on dark matter distribution is displayed in Figure|2] Enhancement to £ 2 induced 
by gas is mild and increases to ~ 1 0% at ~ 0.1/i _1 Mpc, which is consistent with previous works (e.g. 
Jing et all 120061; iGuillet et all l2010h . Gas physics other than adiabatic process does not bring significant 
extra modulation to the two point correlation function of dark matter, their effects are seen in higher order 
functions. 

From Figure |2]it is clear that dark matter clustering is immune to baryons at large scales, stage that gas 
physics play is mainly on sc ales R <~ l/i- x Mpc. In the adiabatic run, at z = 1.442, S° M / S% DM only 
becomes larger than one on scales smaller than l/i^Mpc, it is smaller at z — 0.526, then the boost switches 
to suppression at scales < 0.2/i _1 Mpc. The level of impact is not very large, at 0.1/i _1 Mpc skewness S3 is 
increased by ~ 5% at z = 1.442 and then decreased by ~ 2% at z = 0, variation in kurtosis at such scale 
is ~ +12% at z = 1.442 but 8% at z = 0. 

Effects of star formation activities and other gas physics are much stronger than the adiabatic process, 
amplitudes of 5„s at all redshift are raised significantly, however the increment decreases at lower redshift, 
for instance the relative enhancement to 5*3 and S4 at R = 0.1/i _1 Mpc is ~ 20% and ~ 52% at z = 1.442 
but drops down to ~ 14% and ~ 27% at z — respectively. It appears that if we are about to investigate 
differences among models built with different baryonic processes through clustering analysis of dark matter, 
we should concentrate on higher order statistics on sub-mega parsec scales, and preferentially at higher 
redshift. 



3.2 clustering of baryonic gas 

Correlation functions of gas are presented in Figure |3] Redshift evolution of correlation functions of gas in 
the adiabatic simulation is similar to the dark matter component as in Figure Q] S n being larger at higher 
redshift at scales less than ~ 3/i _1 Mpc. There is more complexity in the non-adiabatic hydrodynamic 
simulation, S n parameters demonstrate intricate evolution path at scales <^ 0.2ft, _1 Mpc, which is probably 
reflection of the composite action from competing gas physical processes, e.g. radiative cooling versus 
feedbacks from supernovae. To separate effect of individual ingred ient of gas physics, seri es simulations 
installed with different prescriptions are definitely needed, like what lvan Daalen et alj d201 ll) did. 

In hydrodynamical simulations, gas distribution deviate from the dark matter distribution in pure dark 
matter more obviously than their dark matter counterpart (Figure|4]i, in aspects of both amplitudes of corre- 
lation functions and affected scale range, while again the effects become we aker at lower red shift. All S n s 
of gas are boosted significantly at scales <^ 4/i~ 1 Mpc, in agreement with Jing et al] (120061) variance of 
gas £ 2 is depressed severely on fairly broad scales extending to around 10/i~ 1 Mpc, at z — £ 2 decreases 
by ~ 30% in adiabatic simulation while by ~ 60% in non-adiabatic simulation at 0.1/i~ 1 Mpc. Bifurcation 
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Fig. 2 Influence of gas physics on dark matter distribution: ratios of correlation functions of dark 
matter in hydrodynamical simulations to those of dark matter in pure dark matter simulation. 



due to differences in gas physical mechanisms employed in simulations is observed too, non-adiabatic gas 
physics inducing stronger variation to the clustering of matter, there is ~ 84% and ~ 220% promotion 
at z = to S3 and S4 respectively at 0.1/i _1 Mpc but only moderate ~ 7% and ~ 14% fortification in 
adiabatic simulation. 



3.3 gas-dark matter clustering segregation 

It is clear now that including gas can result in distribution patterns of both dark matter and gas different to 
the case of pure dark matter at very small scales. From results shown above it appears that gas distribution 
is more severely affected than dark matter, which is easy to understand as all gas physics directly act on gas, 
whilst dark matter is influenced only through its gravitational coupling to the baryons. Direct comparison of 
correlation functions of gas with dark matter is in Figure we can see that for two point correlation function 
differences appear at scales as large as ~ 10/i _1 Mpc but for for higher order cumulants the departure scale 
is ~ 4/i~ 1 Mpc. Volume averaged two point correlation function of gas is obviously lower than the dark 
matter, at 0.1/i _1 Mpc the gap could be as large as ~ 40 — 60%, in adiabatic run the difference at z = is 
around 20% already at l/i _1 Mpc and in the non-adiabatic run it becomes ~ 40% at such scale. However 
S n parameters of gas are much larger than dark matter, for example at 0.1/i _1 Mpc enhancement to the 
skewness S3 at z = is about 10% in adiabatic run and ~ 62% in the non-adiabatic run. The results 
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Fig. 3 £ 2 and S n of gas in hydrodynamic simulations. 



indicate that distribution of gas has a longer tail than that of dark matter, in another words, there are more 
highly concentrated small clumps of gas than dark matter. 

Biasing of gas to dark matter in S„s at z = 0.526 is largest in both hydrodynamic simulations. We 
conjecture that effects of gas physics must reach summit at some redshift z > and then relax after that, 
entering a more passive evolution stage. If galaxy assembly rate is strongly correlated with accumulated 
effect of gas physics, likely there would be an apex within that time interval, of course the exact peak time 
may depends on how baryon processes are cooked. 



4 SUMMARY AND DISCUSSION 



In this paper correlation functions upto the sixth order are estimated from count-in-cell analysis of dark 
matter and gas respectively in three simulations, the pure dark matter run, the run with adiabatic gas process, 
and the one with star formation activities and other gas physics. Major results about influence on matter 
clustering are in the following. 

(1) Compared with the case of pure dark matter, baryon physical processes introduce non-negligible mod- 
ulation to clustering of dark matter, affected regime for dark matter is at scales less than l/i _1 Mpc. 
Adiabatic process alone strengthens £ 2 by ~ 10% at scale 0.1/i _1 Mpc, which is insensitive to redshift; 
S„ parameters in the run deviate from pure dark matter results rather mildly, at 0.1/i _1 Mpc skewness 
S3 evolves from ~ 5% lifting at z = 1.442 to ~ 2% falling at z = 0, meanwhile the difference in 
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Fig. 4 Comparison of correlation functions of gas in hydrodynamic simulations to those of dark 
matter in the pure dark matter run. 



kurtosis S4 changes from ~ +12% to negative ~ —8%. In the run with dissipative gas processes £ 2 
does not differ much to the adiabatic run, but S n parameters all are increased significantly, bringing 
<~ +14% to S3 and ~ +27% to S4 at 0.1ft- _1 Mpc and z = 0, and the amplitude of change is larger in 
higher redshift. 

(2) Gas distribution in hydrodynamic simulations is much more strongly modified than dark matter com- 
ponent. Two point correlation function of gas at z — decreases by ~ 30% in adiabatic simulation 
while by ~ 60% in non-adiabatic simulation at 0.1/i _1 Mpc, the attenuation is weaker at larger scales 
but still obvious at ~ 10/i _1 Mpc. S n parameters of gas are biased upward at scales <^ 4/i~ 1 Mpc, 
dissipative processes add prominently more power to them, giving a ~ 84% promotion at z = to S3 
at 0.1/i _1 Mpc against the moderate ~ 7% fortification in adiabatic simulation. 

(3) There is clustering segregation between gas and dark matter in the same simulation. £ 2 of gas is already 
lower than dark matter counterpart at ~ 10/i _1 Mpc, which is down at 0.1/i _1 Mpc by ~ 40% and 
<~ 62% in the adiabatic run and the non-adiabatic run respectively. S n s of gas are larger than dark 
matter at scales < 4/i _1 Mpc, S3 of gas in adiabatic run is leveled up by ~ 10% while by ~ 60% 
in non-adiabatic run at 0.1/i _1 Mpc. Biasing of gas to dark matter is much stronger in non-adiabatic 
simulation than the adiabatic only run, and the maximal bias is achieved at certain redshift z > 0. 
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Fig. 5 Differences between distributions of gas and dark matter in hydrodynamic simulations. 



It is shown in this work that difference in distribution of dark matter originated from various mecha- 
nisms of gas physics can not be effectively distinguished at the second order level, though apparent discrep- 
ancy appear in gas. It would benefit those applications which rely on second order statistical properties of 
dark matter only, but once going to higher orders one has to consider the systematics brought forward by 
gas. 

Biasing of gas to dark matter is a more interesting problem, aside from that it may be a serious chal- 
lenge to precision c osmology such as the modeling to Sunyaev-Zeldovich effects ("e.g. IShaw et all [20 10; 
lBattagliaetaill2010h . We know that galaxies are biased tracers of dark matter distribution, but galaxies are 
in fact products of gas physics, probably it is more reasonable to assume that galaxies are actually tracing 
gas instead of the dark matter. We conjecture that by the decomposition stochasticity and nonlinearity of 
galaxy bias would be greatly reduced. Standard metho ds exploring relation betwee n galaxies and their host 
halos, such as the halo oc cupation distribut ion model dBeriind & W einberg. 2002) and the conditional lu- 
minosity function model dYang et all 12003). generally use the two point correlation function summarized 
from pure dark matter simulation as reference to the measured galaxy two point correlation function. Data 
points of galaxy two point correlation function usually are at scales from ~ 0.1/i _1 Mpc to a few mega 
parsec within which unfortunately the matter distribution underlying to galaxies is not the same as what 
is in pure dark matter universe, we might have to quantify the amplitude of this kind of systematical bias 
before presenting estimation of number of a particular type of galaxies in halos. 
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